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ABSTRACT  We  report  the  thermal  conductance  G  of  Au/Ti/graphene/Si02  interfaces  (graphene  layers  1  <  n  <  1 0)  typical  of  graphene 
transistor  contacts.  We  find  G  &  25  MW  rrr2  K_1  at  room  temperature,  four  times  smaller  than  the  thermal  conductance  of  a  Au/Ti/ 
Si02  interface,  even  when  n  =  1 .  We  attribute  this  reduction  to  the  thermal  resistance  of  Au/Ti/graphene  and  graphene/Si02  interfaces 
acting  in  series.  The  temperature  dependence  of  G  from  50  <  T  <  500  K  also  indicates  that  heat  is  predominantly  carried  by  phonons 
through  these  interfaces.  Our  findings  suggest  that  metal  contacts  can  limit  not  only  electrical  transport  but  also  thermal  dissipation 
from  submicrometer  graphene  devices. 

KEYWORDS  Interfacial  thermal  conductance,  few-layer  graphene,  metal  contacts,  acoustic  phonons,  heat  dissipation,  thermal 
management  of  graphene  devices 


Graphene,  a  monolayer  or  few  layers  of  sp2-bonded 
carbon  atoms,  has  attracted  immense  attention  over 
the  past  few  years  partly  due  to  the  prospect  of 
replacing  III— V  semiconductors  as  the  workhorse  of  high- 
frequency  electronics.1  In  addition  to  superior  electrical 
properties  (mobility  ~  10000  cm2  V-1  s-1  at  room  temper¬ 
ature2  on  Si02  at  practical  carrier  concentrations  ~1012 
cm-2),  graphene  also  has  exceptionally  high  in-plane  thermal 
conductivity3,4  on  the  order  ~1000  W  m-1  K_1,  which  is 
advantageous  for  thermal  management.  High  in-plane  ther¬ 
mal  conductivity,  however,  implies  that  heat  dissipation 
from  graphenes  is  also  affected  by  how  heat  flows  across 
its  interfaces. 

In  this  Letter,  we  report  the  thermal  conductance  G  of  Au / 
Ti/graphene/Si02  interfaces  which  are  typically  found  in 
graphene  devices,5  for  number  of  graphene  layers  1  <  n  < 
10  and  a  temperature  range  of  50  <  T  <  500  K.  We  find  G 
rs  25  MW  m“2  K-1  irrespective  of  n  at  room  temperature 
and  that  heat  flow  across  the  metal/graphene/Si02  interface 
is  limited  by  G  of  the  metal/graphene  interface  rather  than 
the  graphene/Si02  interface.  Thus,  the  choice  of  metal 
contacts  affects  both  electrical6  and  thermal  transport5  in 
graphene  devices. 

While  in-plane  heat  transport  in  graphene  can  be  at¬ 
tributed  to  acoustic  vibrational  modes  called  phonons,4,7  the 
dominant  heat  carriers  across  graphene  interfaces  in  the 
cross-plane  direction  are  not  known.  Even  under  the  absence 
of  an  external  electric  field,  graphene  could  have  substantial 
carrier  density6  on  the  order  of  1 01 3  cm-2  due  to  charge  transfer 
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from  metal  contacts.8  Since  the  density  of  free  carriers  (elec¬ 
trons  or  holes)  is  high  on  both  sides  of  the  interface,  heat 
transport  across  a  metal/graphene  interface  could  be  facilitated 
by  free  carriers  instead  of  phonons,  analogous  to  heat  transport 
by  electrons  across  metal/metal  interfaces.9  Moreover,  since 
carriers  in  graphene  are  readily  scattered  by  polar  optical 
phonons  (sometimes  called  remote  interfacial  phonons  or 
surface-optical  phonons)  in  the  Si02  substrate,2,10  heat  trans¬ 
port  through  graphene/Si02  interfaces  could  have  a  contribu¬ 
tion  from  direct  energy  transfer11  between  free  carriers  in 
graphene  and  polar  optical  phonons  in  Si02. 

To  identify  the  dominant  heat  carriers  across  graphene/ 
metal  and  graphene/Si02  interfaces,  we  measured  the  ther¬ 
mal  conductance  G  of  Au/Ti/graphene/Si02  interfaces  over 
a  temperature  range  of  50  <  T  <  500  K.  We  find  that  G  of 
such  interfaces  depends  only  weakly  on  temperature  above 
1 00  K.  The  temperature  dependence  of  G  is  expected  to 
follow  the  temperature  dependence  of  heat  capacity  for  the 
excitations  that  are  most  responsible  for  heat  conduction.12 
The  heat  capacity  of  the  degenerate  electrons  in  Au  and  the 
heat  capacity  of  optical  phonons  in  Si02  increase  signifi¬ 
cantly  between  1 00  and  500  K;  the  contribution  of  acoustic 
phonons  to  the  heat  capacity  of  Au,  on  the  other  hand,  is 
essentially  constant  in  this  temperature  range  due  to  the  low 
Debye  temperature  of  Au.  Therefore,  we  conclude  that 
acoustic  phonons  are  the  dominant  carriers  of  heat  across 
these  interfaces. 

Even  when  an  interface  is  atomically  abrupt,  heat  trans¬ 
port  by  phonons  across  it  is  limited  by  the  finite  probability 
of  phonon  transmission;  the  resulting  impedance  to  heat 
flow  is  characterized  by  a  property  of  an  individual  interface 
called  the  Kapitza  thermal  conductance,12  Gh  The  total 
thermal  conductance  per  unit  area  G  of  a  thin  film  sand- 
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wiched  between  two  solids  can  often  be  attributed  to  the 
series  sum  of  contribution  from  two  discrete  interfaces  and 
the  bulk  of  the  thin  film.13  If  this  assumption  were  valid  for 
the  case  of  metal/graphene/Si02  interfaces,  we  could  write 

G  )metal/n-LG  ^1  )n-LG/Si02  +  (/i/A)n.LG  (1) 


where  h  and  A  are  the  thickness  and  the  cross-plane  thermal 
conductivity  of  n-layer  graphene  (n-LG),  respectively,  and 
hi  A  is  the  corresponding  thermal  resistance  per  unit  area. 
Assuming  A  >  1  W  mT1  K_1,  hi  A  <  10“9  m2  K  and  is 
negligible  compared  to  Gf  b 

When  the  him  is  sufficiently  thin,  however,  phonons 
(especially  long-wavelength  phonons)  could,  in  principle, 
traverse  directly  between  the  metal  and  SiC>2  without  being 
scattered  by  the  interfaces.14  In  such  cases,  eq  1  is  invalid, 
and  we  should  treat  the  structure  as  one  diffuse  interface 
instead  of  two  discrete  interfaces.  Factors  dictating  at  what 
him  thickness  eq  I  breaks  down  are  not  yet  established. 
However,  from  our  measurements  on  G  of  Au/Ti/graphene/ 
Si02  interfaces,  we  hnd  that  G  remains  unchanged  for 
number  of  graphene  layers  of  1  <  n  <  10  (equivalent  to 
thickness  of  0.4  <  h  <  4  nm).  Moreover,  the  measured  G  is 
comparable  to  the  G  calculated  using  eq  1  over  a  tempera¬ 
ture  range  of  50  <  T  <  500  K.  Thus,  our  results  suggest  that, 
even  when  the  graphene  layer  is  only  one  atomic  layer  thick, 
heat  transport  across  the  metal/graphene/oxide  interfaces 
is  governed  by  the  thermal  conductance  of  two  (metal/ 
graphene  and  graphene/oxide)  uncoupled  interfaces.  This 
indicates  that  heat  conduction  due  to  direct  transmission  of 
phonons  between  Au  and  Si02  is  negligible. 

Our  test  structures  resemble  typical  metal  contacts  in 
graphene  devices;5  see  Figure  la.  Prior  to  deposition  of 
graphene,  we  cleaned  selected  Si02/Si  substrates  by  1 0  min 
of  sonication  in  acetone  and/or  heat  treatment  at  200  °C  in 
ambient  for  1  h;  we  noticed  no  difference  in  the  measured 
G  with  or  without  precleaning.  We  deposited  graphene  on 
104  nm  Si02  on  Si  by  micromechanical  exfoliation15  from 
natural  graphite.  We  then  annealed  the  substrate  in  a 
chemical  vapor  deposition  (CVD)  furnace  at  400  °C  for  35 
min  in  Ar/H2  mixture  gas  to  remove  the  adhesive  residue 
from  the  tape.  We  patterned  «100  nm  thick  Au  pads  on 
selected  graphene  flakes  by  e-beam  lithography,  e-beam 
evaporation,  and  lift-off  with  ~2  nm  Ti  predeposited  as  the 
adhesion  layer.  E-beam  evaporation  was  performed  at  a 
base-pressure  of  7  x  1CT7  Torr,  and  the  Au  and  Ti  were 
deposited  at  rates  of  2.5  and  0.5  A/s,  respectively.  Due  to 
low  deposition  rate  of  Ti,  a  significant  amount  of  oxygen 
could  be  trapped  in  the  Ti  film.  We  note  that  even  though 
the  metal  films  were  deposited  in  vacuum,  water  vapor  that 
adsorbs  on  the  substrate  during  the  air  exposure  after  the 
annealing  step  may  still  be  present  during  the  deposition  of 
the  Au/Ti  contacts.  We  purposely  deposited  large  metal  pads 
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FIGURE  1 .  (a)  Vertical  cross  section  of  the  samples  (not-to-scale). 
Au(100  nm)/Ti(2  nm)  thin  metal  pads  were  patterned  on  graphene 
flakes  on  Si02  (104  nm)  on  Si  substrate  by  e-beam  lithography,  (b) 
Example  optical  microscopy  image  of  a  graphene  sample  covered 
by  metal  pads. 


on  our  samples  such  that  the  same  metal  pad  covers  both 
regions  with  and  without  graphene  flakes;  see  Figure  lb. 
Time-domain  thermoreflectance  (TDTR)  measurements,  dis¬ 
cussed  below,  were  performed  on  both  regions  with  and 
without  graphene  flakes  covered  by  the  same  metal  pads  to 
derive  G  of  Au/Ti/graphene/Si02  interfaces. 

Our  samples  consist  of  graphenes  with  number  of  layers 
from  1  <  n  <  10,  denoted  by  n-LG  (n  layers  of  graphene) 
throughout  this  work.  We  employed  a  newly  developed 
approach  based  on  Raman  spectroscopy16  to  count  the 
number  of  layers  n  of  the  graphene  flakes.  In  this  approach, 
n  is  determined  from  the  ratio  of  the  integrated  intensity  of 
Raman  peak  at  ~520  crrr1  due  to  scattering  by  optical 
phonons  in  silicon,  /(Si),  and  integrated  intensity  of  G  peak 
due  to  scattering  by  doubly  degenerate  in-plane  optical 
phonons  (iTO  and  LO)  in  graphene,  /(G);  see  Figure  2a.  We 
validated  n  derived  using  this  approach  by  measuring  the 
thickness  of  the  graphenes  by  atomic  force  microscopy 
(AFM);  a  more  complete  description  of  this  approach  will  be 
published  in  ref  16.  We  measure  the  Raman  spectra  of 
graphenes  by  a  custom-built  Raman  spectrometer;  a  488  nm 
laser  beam  excites  the  Raman  spectra.  We  use  laser  powers 
of  mW  in  our  measurements  and  a  20 x  objective  lens 
with  N.A.  =  0.4  to  focus  the  laser  beam  and  collect  Raman- 
scattered  light  in  all  polarizations.  The  full  width  half¬ 
maximum  spectral  resolution  of  our  Raman  setup  is 
cirT1 . 

To  verify  that  the  graphene  flakes  are  not  damaged  by 
deposition  of  Au/Ti,  we  coated  a  monolayer  graphene  (1- 
LG)  sample  with  a  semitransparent  layer  of  Au  (8  nm)/Ti  (2 
nm)  and  compared  the  Raman  spectrum  of  the  sample 
before  and  after  metal  deposition;  see  Figure  2b.  We  observe 
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FIGURE  2.  (a)  Determining  graphene  layer  number  (n)  from  ratios 
of  integrated  intensity  of  the  Raman  peak  due  to  scattering  by 
optical  phonons  in  silicon,  /(Si),  and  integrated  intensity  of  graphene 
G  peak,  /(G).  The  Raman  spectra  were  acquired  using  a  488  nm  laser. 
Raman  measurements  (open  circles)  are  compared  to  average  /(G)/ 
/(Si)  values  from  ref  16  (dashed  line)  to  determine  the  number  of 
layers  n.  The  thickness  of  two  selected  graphene  flakes  is  measured 
by  atomic  force  microscopy  (AFM),  as  labeled,  (b)  Raman  spectra  of 
monolayer  graphene  before  and  after  coverage  with  Au  (8  nm)/Ti 
(2  nm).  The  spectrum  of  graphene  after  metal  deposition  is  vertically 
down-shifted  by  170  counts  mW s  1  for  clarity.  D  peaks  (defect- 
mediated  Raman  scattering  from  iTO  phonons  near  the  Brillouin 
zone  boundary)  are  weak  in  both  spectra,  indicating  that  metal 
deposition  does  not  significantly  damage  the  graphene  structure. 
Metal  deposition  causes  the  G  peak  to  split  into  two  slightly 
broadened  peaks  centered  at  1590  and  1622  cm"1. 

a  distinct  2D  peak  (due  to  two  iTO  phonons  near  Brillouin 
zone  boundary)  at  2707  crrT1  with  a  fwhm  of  35  cm-1  after 
the  deposition  of  thin  metal  layers.  This  narrow  2D  peak  at 
~2700  crrT1  is  a  fingerprint17  of  monolayer  graphene  and 
indicates  that  graphene  was  not  significantly  damaged  by 
e-beam  deposition  of  the  metal  film.  The  D  peak  (due  to  one 
iTO  phonon  near  Brillouin  zone  boundary  and  defects)  at 
~  1 350  crrT1  is  weak  after  the  metal  deposition,  suggesting 
that  the  density  of  graphene  defects  remains  low.  The  G  peak 
is  slightly  split  and  broadened  after  metal  evaporation;  see 
Figure  2b.  This  change  of  shape  is  inconsistent  with  a 
sharpened,  red-shifted  G  peak  due  to  Landau  damping  and 
stiffening  of  phonons  by  high  carrier  concentrations  in 
graphene;18  we  are  unsure  of  the  cause  of  these  changes  in 
the  shape  of  the  G  peak  but  interactions  between  the  G-band 
optical  phonons  in  graphene  and  the  electronic  excitations 
of  the  metal  may  be  playing  a  role. 

We  measure  the  thermal  conductance  G  by  time-domain 
thermoreflectance19'20  (TDTR).  A  schematic  diagram  of  our 


equipment  is  given  in  ref  20  and  our  method  for  data 
analysis  is  described  in  ref  2 1 .  In  our  TDTR  measurements, 
a  laser  beam  from  a  mode-locked  Ti:sapphire  laser  is  split 
into  a  pump  beam  and  a  probe  beam  with  the  relative  optical 
path  being  adjusted  via  a  mechanical  delay  stage.  The  pump 
beam  is  modulated  at  a  frequency/=  1 0  MHz.  Laser  pulses 
of  the  pump  beam  of  <1  nj  are  absorbed  by  the  Au  pads, 
inducing  a  temperature  oscillation  of  <10  K  at  the  surface 
of  the  Au  pads.  Cooling  of  the  surface  after  being  heated  by 
pump  pulses  is  then  monitored  through  changes  in  the 
intensity  of  the  reflected  probe  beam  using  a  Si  photodiode 
and  a  lock-in  amplifier.  The  radii  of  the  laser  beams  used 
during  the  measurements  are  either  7.5  or  3.8  fim  at  the 
sample  surface.  We  used  total  laser  powers  of  <11 0  mW  for 
laser  spot  size  of  7.5  fim  and  <45  mW  for  laser  spot  size  of 
3.8  fim,  creating  temperature  rises  of  <10  K.  The  ratios  of 
the  power  of  the  pump  beam  to  the  probe  beam  vary  from 
1:1  to  20;  1 ;  as  expected,  we  did  not  observe  any  difference 
in  the  derived  G  using  different  power  ratios. 

To  enhance  the  accuracy  of  G  of  Au/Ti/n-LG/SiC>2  inter¬ 
faces  derived  using  TDTR,  we  performed  TDTR  measure¬ 
ments  on  adjacent  regions  with  and  without  graphene  flakes 
for  every  graphene  flake  we  studied.  We  compared  the  ratio 
of  in-phase  and  out-of-phase  of  the  lock-in  amplifier  of  both 
measurements  to  numerical  solutions  of  a  thermal  model21 
that  takes  into  account  changing  of  the  radius  of  the  pump 
beam  at  different  relative  delay  times;22  see  Figure  3a.  We 
first  take  G  =  1 00  MW  m"2  K"1  for  Au/Ti/SiC>2  interfaces  (this 
value  of  G  was  prior  measured  using  a  500  nm  Si02  coated 
with  Au/Ti)  and  lump  the  uncertainty  of  TDTR  measure¬ 
ments  by  deriving  the  thermal  conductivity  of  Si02  thin  films 
from  initial  measurements  on  regions  without  graphene 
flakes.  Since  the  uncertainty  should  be  similar  for  measure¬ 
ments  on  adjacent  regions  with  graphene  flakes,  we  then 
derived  G  of  Au/Ti/n-LG/Si02  interfaces  by  fitting  the  second 
TDTR  measurements  on  regions  with  graphene  flakes  using 
the  thermal  conductivity  of  Si02  previously  derived.  The 
uncertainty  of  G  for  Au/Ti/n-LG/Si02  interfaces  derived  using 
this  approach  is  ^20%.  The  main  contributor  to  this  is  the 
uncertainty  (~35  % )  in  G  of  the  Au/Ti/Si02  interfaces,  which 
contributes  »12%  to  the  total  uncertainty. 

We  plot  G  of  Au/Ti/n-LG/Si02  interfaces  as  a  function  of 
number  of  graphene  layer  n  in  Figure  3b.  We  find  that  at 
room  temperature  20  <  G  <  30  MW  m-2  K_1  for  Au/Ti/n-LG/ 
Si02  interfaces,  irrespective  of  n;  this  value  of  G  is  a  factor  of  4 
smaller  than  G  =  100  MW  m“2  K_1  directly  measured  on  our 
control  Au/Ti/Si02  interfaces;  see  Figure  3b.  Similar  values  of 
thermal  conductance  were  reported  for  Au/SAM/GaAs  inter¬ 
faces,23  in  which  a  self-assembled  monolayer  (SAM)  of  al- 
kanedithiols  is  sandwiched  between  Au  and  GaAs. 

We  note  that  this  relatively  high  value  of  G  for  Au/Ti/n- 
LG/Si02  interfaces  implies  that  graphene  capped  with  metals 
cannot  be  partially  suspended  between  asperities  of  Si02 
substrate,  which  is  generally  perceived  to  be  the  case  for 
uncapped  graphenes  supported  on  Si02. 24,25  If  the  capped 
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FIGURE  3.  (a)  Ratios  of  in-phase  and  out-of-phase  signals  (open  circles) 
measured  by  the  lock-in  amplifier  as  a  function  of  delay  time  t  between 
pump  and  probe  pulses.  Measurements  on  Au/Ti/Si02/Si  (no  graphene), 
Au/Ti/ 1  -LG/Si02/Si,  and  Au/Ti/3-LG/Si02/Si  are  labeled  as  “no  graphene”, 
“1-LG”,  and  “3-LG”,  respectively.  The  monolayer  (1-LG)  and  trilayer 
(3-LG)  graphenes  are  deposited  on  the  same  substrate  and  covered  with 
metal  pads  under  the  same  conditions,  as  shown  in  Figure  lb.  The  solid 
lines  are  calculations  of  a  thermal  model,  (b)  Thermal  conductance  G 
per  unit  area  of  Au/Ti/n-LG/Si02  measured  by  TDTR  using  laser  beams 
with  1/e2  spot  size  of  7.5  fim  (open  circles)  and  3.8  f/m  (solid  circles)  as 
a  function  of  number  of  layers  n;  n  =  0  corresponds  to  G  of  the  Au/Ti/ 
Si02  interface.  The  solid  line  is  a  calculation  using  eq  1 .  The  uncertainty 
of  measurements  is  ~20% . 

graphene  were  suspended  between  asperities  of  Si02,  the 
contact  area  between  graphene  and  SiC>2  would  be  small  and 
the  thermal  conductance  of  graphene/Si02  interfaces  ex¬ 
tremely  low.  For  example,  if  graphenes  in  our  measure¬ 
ments  were  partially  suspended  across  voids  with  a  contact 
area  of  <10%,  our  measurements  indicate  that  the  real 
thermal  conductance  of  the  metal/graphene/oxide  is  >250 
MW  rrT2  K~\  inconsistent  with  our  measurement  of  G  of 
an  Au/Ti/Si02  interface  and  prior  measurements  of  G  «  85 
MW  rrT2  K~'  for  the  1  -LG/Si02  interfaces26  and  G  <  1 20  MW 
rrT2  FT'  for  the  metal/graphite  interfaces.27  Thus,  our  results 
indicate  that  either  the  capped  graphene  conforms28  to  the 
Si02  substrate  or  the  voids  are  filled  with  water  or  hydro¬ 
carbons  that  can  carry  heat  across  the  gap. 

Our  assertion  that  the  encased  graphenes  conform  to  the 
substrate  is  reinforced  by  our  time-domain  pump-beam- 
deflection29  (TD-PBD)  measurements  of  the  acoustic  waves 
reflected  by  interfaces  with  and  without  graphenes.  We  find 
that  the  shape  and  amplitude  of  the  acoustic  reflections  are 
comparable  with  and  without  graphenes.  Since  voids  with 
zero  acoustic  impedance  should  vastly  enhance  the  acoustic 
reflection,  our  findings  suggest  that  the  area  occupied  by 
voids  is  a  small  fraction  of  the  total  area. 
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Our  values  for  G  of  graphene  contacts  suggest  that  heat 
dissipation  from  graphene  field-effect  transistors  (GFETs) 
with  dimensions  I  <  600  nm  will  often  be  affected  by  the 
thermal  conductance  of  graphene  interfaces.  To  estimate  the 
importance  of  graphene  interfaces  in  the  thermal  manage¬ 
ment  of  GFETs  supported  on  Si02/Si,  we  assume  that  heat 
is  generated  at  the  center  of  the  GFET  and  can  dissipate  to 
the  surroundings  either  through  metal  contacts  or  the  Si02 
substrate.  The  dominant  path  of  heat  dissipation  is  governed 
by  a  characteristic  length  called  the  thermal  healing  length30,31 

Zh  =  VhA/G  (2) 

where  h  =  0.35  nm  and  A  as  1000  W  rrT1  K_1  are  the 
thickness  and  in-plane  thermal  conductivity  of  monolayer 
graphene  and  1/G  =  1  /Gi.LG/sio2  +  hsioJ^sio2-  For  monolayer 
graphene  supported  on  hS\o2  =  100  and  300  nm,  Zh  aa  180 
and  290  nm,  respectively.  For  L  <  2 Zh,  heat  dissipates  from 
the  graphene  mainly  through  the  metal  contacts;  while  for 
L  >  2/h,  heat  dissipates  mainly  through  the  Si02  substrate. 
Although  most  current  graphene  devices  have  I  >2/h,  the 
characteristic  dimensions  of  future  graphene  devices  will  be 
reduced  <2lh  leading  to  heat  dissipation  mainly  through  the 
metal  contacts. 

To  estimate  the  importance  of  G  compared  to  other 
thermal  resistances  of  a  GFET,  we  consider  a  hypothetical 
device  of  length  I  =  100  nm  between  contacts,  width  W  = 
400  nm,  on  /isl0,  =  100  nm,  and  100  x  400  nm2  metal 
contact  area.  The  estimated  heat  flow  resistances  through 
one  contact,  along  the  graphene,  and  through  the  Si02  are 
Rc^  106  K  W_1,  Rc  «  7.1  x  105KW“',  and  Rsio2^2.2  x 
1 06  K  W_1 ,  respectively.  While  heat  dissipation  from  GFETs 
will  depend  on  the  realistic  three-dimensional  geometry  of 
the  device,  this  simple  estimate  reveals  the  heat  flow  path 
through  the  metal  contacts  could  in  fact  be  significant  for 
thermal  management  of  future  graphene  devices. 

To  further  enhance  our  understanding  of  heat  conduction 
across  graphene,  we  measured  G  of  Au/Ti/n-LG/Si02  inter¬ 
faces  over  a  wide  temperature  range  of  50  <  T  <  500  K; 
see  Figure  4.  As  introduced  above,  the  lack  of  temperature 
dependence  of  G  for  T  >  1 00  K  as  in  Figure  4  supports  our 
assertion  that  the  dominant  heat  carriers  across  Au/Ti/n-LG 
and  n-LG/Si02  interfaces  are  acoustic  phonons  and  not  free 
carriers9  or  direct  coupling1 1  of  free  carriers  in  graphene  and 
optical  phonons  in  Si02.  Our  conclusion  that  free  carriers 
are  not  responsible  for  heat  conduction  across  graphene 
interfaces  is  also  consistent  with  an  estimate  using  the 
interface  form  of  the  Wiedemann— Franz  law.9  For  agraphene/ 
metal  contact  area  A  =  1 00  /im2  with  an  electrical  resis¬ 
tance32  Re  as  10  £2,  the  thermal  conductance  at  room 
temperature  is  Ge  «a  LTI(ARe)  10  kW  rrT2  K_1  (Lorenz 
number  L  =  2.45  x  1 0~8  £2  W  K-2),  a  factor  of  1 000  lower 
than  the  measured  lattice  thermal  conductance  of  >1 0  MW 
rrT2  K~'.  In  other  words,  except  for  superb  electrical  contacts 
that  are  yet  to  be  discovered,  heat  transport  across  metal/ 
graphene  is  predominantly  carried  by  phonons. 
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FIGURE  4 .  Temperature  dependence  of  thermal  conductance  G  of  Au/ 
Ti/1-LG/Si02  (solid  triangles)  and  Au/Ti/3-LG/Si02  (solid  circles)  inter¬ 
faces.  Thermal  conductance  of  Au/Ti/Si02  (solid  diamond,  this  work), 
Au/Ti/graphite  (solid  square,  this  work),  1-LG/Si02  from  Chen  et  al. 
(open  right  triangles,  ref  26)  and  Freitag  et  al.  (open  triangle  down,  ref 
11),  Au/l-LG  (open  triange  up,  Cai  et  al.,  ref  37),  SWCNT/Si02  (open 
circle,  Pop  et  al.,  ref  31),  Au/graphite  (open  squares,  ref  27),  and  Ti/ 
graphite  (open  diamonds,  ref  27)  are  included  for  comparison.  SWCNT 
denotes  single-wall  carbon  nanotube.  The  dashed  lines  are  diffuse 
mismatch  model12  (DMM)  calculations  for  G  of  Au/Ti/Si02  (top)  and  Au/ 
Ti/graphite  (bottom)  interfaces.  The  solid  line  is  the  prediction  of  eq  1 
using  DMM  calculations  of  GAu/Ti/graphite  and  measurements  of 
Gi-LG/sio2  from  ref  26.  Equation  1  agrees  with  measurements  of  Au/Ti/ 
n-LG/Si02  over  the  entire  temperature  range. 

The  reduction  of  thermal  conductance  due  to  the  pres¬ 
ence  of  atomic  graphitic  layers  that  we  observe  in  Figure  3b 
could  be  explained  by  weaker  coupling33  between  metal  and 
oxide.  Under  this  explanation,  the  metal/graphene/oxide 
structure  is  treated  as  a  perturbation  that  alters  the  prob¬ 
ability  for  phonons  to  transmit  between  metal  and  oxide, 
instead  of  being  treated  as  two  discrete  interfaces.  We  argue, 
however,  that  this  explanation  is  inadequate  especially  for 
thick  graphenes  of  n  =  10  with  layer  thickness  of  ~4  nm. 
Since  most  heat  is  carried  by  phonons  with  wavelength  <1 
nm  at  room  temperature,  a  4  nm  layer  is  too  thick  to  be 
considered  a  diffuse  interface.  Moreover,  most  phonon 
modes  in  metal  and  oxide  have  relatively  low  energy  (the 
highest  energy  of  acoustic  phonons  is  1 8  and  62  meV  in  Au 
and  amorphous  Si02,  respectively),  but  the  density  of  low- 
energy  modes  in  graphene  is  small.  To  give  an  idea  of  the 
density  of  low-energy  modes  available  in  these  materials, 
the  ratio  of  heat  capacity  of  Au,  graphite,  and  Si02  at  30  K  is 
22:1:2.6.  Consequently,  phonons  cannot  traverse  directly 
and  coherently  between  metal  and  oxide. 

Instead,  we  attribute  the  lower  thermal  conductance  of 
Au/Ti/n-LG/Si02  interfaces  to  Kapitza  thermal  resistance12  of 
Au/Ti/n-LG  and  n-LG/Si02  interfaces  acting  in  series.  To  test 
this  hypothesis,  we  compare  our  measurements  to  G  calcu¬ 
lated  using  eq  1 ;  see  the  solid  line  in  Figure  3b.  In  the 
calculation,  we  use  a  prior  measurement26  of  G  of  a  1  -LG/ 
Si02  interface  for  Gi-Lg/sio2  and  approximate  Gmetai/i-LG  from 
our  measurement  of  G  of  Au/Ti/graphite  interfaces.  We  find 
that  the  calculation  of  eq  1  is  in  reasonable  agreement  with 
our  measurements  of  G  of  Au/Ti/n-LG/Si02  interfaces.  Our 
finding  suggests  that  even  when  graphene  is  only  one  atomic 


layer  thick,  heat  is  carried  across  the  metal/graphene/oxide 
interfaces  by  two-stage  transmission  of  acoustic  phonons 
across  metal/graphene  and  graphene/oxide  interfaces. 

To  further  support  our  assertion  that  heat  transport  across 
the  interface  is  dominated  by  the  Kapitza  thermal  resistance 
of  two  graphene  interfaces,  we  compare  our  G  measure¬ 
ments  over  a  wide  temperature  range  to  G  calculated  using 
eq  1  in  Figure  4.  In  the  calculations,  we  use  prior  measure¬ 
ments26  of  a  1  -LG/Si02  interface  for  Gi.uysioz  over  the  similar 
temperature  range.  Since  the  temperature  dependence  of 
G  of  Au/Ti/graphene  is  not  available,  we  approximate 
Gmetai/i-LG  from  G  of  Au/Ti/graphite  interfaces  calculated  using 
the  diffuse  mismatch  model  (DMM).12  Our  implementation 
of  the  DMM  calculations  is  slightly  different  from  the  ap¬ 
proach  described  in  ref  12,  which  is  derived  for  G  at  low 
temperatures.  Following  a  previous  approach  to  calculations 
of  thermal  conductivity22  at  high  temperatures,  we  assume 
a  linear  dispersion  for  phonons  in  Au  and  set  cutoff  frequen¬ 
cies  using  the  frequencies  of  longitudinal  and  transverse 
acoustic  phonons  at  the  Brillouin  zone  boundary  in  the  [1 00] 
direction;  the  cutoff  frequencies  for  longitudinal  and  trans¬ 
verse  phonons  in  Au  are  4.4  and  2.8  THz,  respectively.  We 
assume  that  scattering  at  the  interfaces  is  elastic  and  allow 
mode  conversion  at  the  interfaces  as  long  as  the  energy  is 
conserved.  We  follow  ref  12  and  approximate  the  transmis¬ 
sion  coefficient  a  of  phonons  from  Au  to  Si02  using 

4io2  ^vj  )sio2 

a  =  7 - =  - =7 - (3) 

si°2  au  (Sia  )si0^  +  (Evj  )Au 

where  /Au  =  (2Vj~2)Au  is  the  sum  of  vf2  of  all  phonon  modes 
j  available  in  the  Au  side  and  Vj  is  the  speed  of  sound  of 
phonons  with  mode  j.  We  note  that  IAu  depends  on  the 
phonon  frequency  v  due  to  the  cutoff  frequencies  that  we 
impose.  For  Au  to  Si02,  a  =  0. 1 0  when  v  <  2.8  THz  and  a  = 
0.65  when  v  >  2.8  THz.  However,  the  transmission  coef¬ 
ficient  of  phonons  from  Au  to  graphite  could  not  be  esti¬ 
mated  using  eq  3  since  graphite  is  highly  anisotropic.  We 
therefore  derive  the  effective  /graphite  =  6.25  x  1CT8  s2  rrT2 
by  fitting  the  calculations  to  a  measurement  of  G  =  52  MW 
rrr2  for  an  Au/Ti/graphite  interface  at  room  temperature; 
see  Figure  4.  The  derived  a  for  Au  to  graphite  is  then  a  = 
0.04  when  v  <  2.8  THz  and  a  =  0.40  when  v  >  2.8  THz.  Our 
approach  of  deriving  a  from  fitting  is  compatible  to  the 
approach  by  ref  34,  which  estimates  a  from  the  two- 
dimensional  density  of  states  in  graphite. 

Our  measurements  of  G  agree  well  with  the  calculations 
using  eq  1  over  the  entire  temperature  range,  see  Figure  4, 
indicating  that  heat  flow  across  graphene  is  governed  by  the 
Kapitza  thermal  conductance  of  the  two  interfaces  acting  in 
series.  Our  finding  that  the  bulk  of  graphene  does  not 
significantly  contribute  to  total  thermal  resistance  of  graphene 
in  the  cross-plane  direction  is  consistent  with  similar  conclu- 
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sions  reached  for  heat  flow  across  carbon  nanotubes,3 1,35 
SAMs,23  and  molecular  chains,36  which  are  also  limited  by 
interfaces  and  thus  are  insensitive  to  the  size  of  the  nano¬ 
tubes,  SAMs,  or  the  molecular  chains.  Similar  results  are 
observed  in  graphene;  see  Figure  3b. 

In  conclusion,  we  find  that  heat  conduction  across  metal/ 
graphene/oxide  interfaces  is  limited  by  finite  transmission 
of  phonons  between  metal  and  graphene.  Thus,  for  thermal 
management  of  graphene  devices,  metals  with  high  Debye 
temperature  (e.g.,  Cr,  Ni,  Ti,  Al)  could  be  better  choices  of 
metal  thermal  contacts  because  of  better  energy  match 
between  phonon  modes  in  these  metals  and  in  graphene. 
Our  results  enhance  the  microscopic  understanding  on  heat 
transport  across  a  single  layer  of  crystalline  atoms. 
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